Adiabatic Condition for Nonlinear Systems 
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We revisit the adiabatic criterion in stimulated Raman adiabatic passage for the three-level A- 
system, and compare the situation with and without nonlinearity. In linear systems, the adiabatic 
condition is derived with the help of the instantaneous eigenvalues and eigenstates of the Hamil- 
tonian, a procedure that breaks down in the presence of nonlinearity. Using an explicit example 
relevant to photoassociation of atoms into diatomic molecules, we demonstrate that the proper way 
to derive the adiabatic condition for the nonlinear systems is through a linearization procedure. 
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According to the adiabatic theorem of quantum me- 
chanics [lj, if the Hamiltonian is changed sufficiently 
slowly, then a system in a given non-degenerate eigen- 
state of the initial Hamiltonian (say, \i)) evolves into the 
corresponding eigenstate of the instantaneous Hamilto- 
nian without making any transitions. Here "sufficiently 
slowly" means that the rate of change of the Hamiltonian 
is much smaller compared to the level spacings: 



|ff | < h\uo fl \ 2 , (for any / ^ i) 



(1) 



where u)u denotes the transition frequency between the 
instantaneous eigenstates \i(t)} and \f{t)}. Condition JTJ 
is referred to as the adiabatic condition. 

The standard version of the adiabatic theorem, which 
has many applications in quantum state preparation and 
manipulation, applies to linear quantum systems. We, 
however, have noticed that the adiabatic condition |QJ 
have been used in recent studies of nonlinear quantum 
gases. This may be problematic since the proof of the 
adiabatic condition makes explicit use of the concept of 
an orthonormal set of energy eigenstates and the lin- 
ear superposition principle involving these states, both 
of which become invalid when nonlinearity is introduced 
into the system. The purpose of this paper is to pro- 
vide a general method for deriving the adiabatic condi- 
tion in nonlinear systems. This is achieved by consider- 
ing a specific example of coherent population transfer in 
three-level quantum systems using the stimulated Raman 
adiabatic passage (STIRAP) method 0,0. 

Consider the three-level A-system schematically shown 
in Fig.^ The excited state |e) is coupled to two ground 
states | a) and \g), with the coupling strengths (Rabi fre- 
quencies) denoted as Jl p and f^, respectively. Here the 
subscripts "p" ( "cf ' ) stand for "pump" ( "dump" ) . In the 
interaction picture, under the two-photon resonance con- 
dition (<5 = 0), the Hamiltonian of this system reads 

H lin = -hA\e)(e\ + ^ (Q p \a)(e\ + tl d \g)(e\ + h.c.) . (2) 

Without loss of generality, we will assume that the Rabi 
frequencies fl p ,d are real and positive. 
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FIG. 1: The three-level system. Q p and £ld are the two cou- 
pling strengths, A and S are one and two-photon detunings, 
respectively. 



Hamiltonian J2J can be easily diagonalized. The par- 
ticular energy eigenstate 



|CPT) = — |o> - O p \g)) 

"off 



where f2 e fi = y O 2 + f^, is known as the coherent popu- 
lation trapping (CPT) state or the dark state. The state 
|CPT) has zero eigenenergy and is decoupled from the 
excited state. When A = 0, the other two eigenstates 
possess energies ±hQ, e g/2. When Q, Pj 4 are varied adia- 
batically, an initially prepared CPT state will remain in 
the instantaneous CPT state. A straightforward appli- 
cation of 0) leads to the adiabatic condition [2j 



Hin = 



0: 



\X\ 



- L «l, (3) 



off 



where x = Op/Od- This property of the CPT state fa- 
cilitates a coherent population transfer from \a) to \g) 
when x changes from to oo. This is the basis for STI- 
RAP which has wide applications ranging from chemical- 
reaction dynamics Q to laser cooling of atoms . 

In this paper, we want to derive the correct adiabatic 
condition for a nonlinear three-level system describing a 
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coupled atom-molecule system. This system has recently 
received significant attention due to the possibility of cre- 
ating ultracold molecules by associating cold atoms 0,0- 
For the nonlinear system, we use the same level scheme 
of Fig. ^ Now state \a) represents an atomic state, while 
|e) and \g) are excited and ground diatomic molecular 
states. To keep things as simple as possible without sac- 
rificing essential physics, we neglect nonlinear collisions 
between particles. The only nonlinearity comes from the 
fact that it takes two atoms to form a molecule. Under 
the two-photon resonance condition, the Hamiltonian in 
second quantized form reads 



fiA^l^e + ^ (flplpllpalpa + fitful 



h.c. 



(4) 

where tpi and tpj are the annihilation and creation oper- 
ators for state \i), respectively. In the mean-field treat- 
ment, these are replaced by c- numbers tpi and ip* , which 
obey the following coupled equations: 



ilpa = £l p 1p* a 1pe , 

ilpe = Atp e + ^-Ipl + ^tpg , 
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fi* 



-Ipe 



(5a) 
(5b) 

(5c) 

tpi 8 are normalized as \ip a \ 2 + 2(|V'g| 2 + IV'el 2 ) = 1, a 
consequence of the conservation of total atom numbers. 

To find the eigenstates of the system, we simply replace 
id/dt at the left hand sides of Eqs. © by u>, the eigenfre- 
quency. It is easy to show, as in the linear counterpart, 
the nonlinear A-system supports a CPT eigenstate with 
zero eigenenergy 7] , with the corresponding state vector 
given by = (ip°, tp®, tpg) T where 



2fi d 



fi„ 



^eff 



1/2 



(6) 



with nf s 



fi 2 



8^. 



Evidence of the CPT state 

in coupled atom-molecule systems has been reported in 
several recent experiments |8j . 

Focusing on the case of vanishing single-photon detun- 
ing, A = 0, we can easily obtain the other eigenstates: 

^(0,±l,lf, 

with eigenfrequencies ifi^/2, respectively. Furthermore, 
when fi^/fip < 1, two more eigenstates appear: 




with eigenfrequencies ±fi p /2, respectively. That there 
may exist more eigenstates than the dimension of the 
Hilbert space is unique for nonlinear systems Q. The 
CPT state distinguishes itself from others with vanishing 



probability amplitude in the excited state ip e . The non- 
orthorgonality between the CPT state and any of the 
other eigenstates is quite obvious. 

With the lack of a set of orthonormal energy eigen- 
states, the conventional linear adiabatic condition Q or 
© can no longer be used for the nonlinear system. We 
want to derive the right adiabatic condition suitable for 
the nonlinear A-system under study. Our goal is to find 
the condition under which the system stays in the instan- 
taneous CPT state © when the Rabi frequencies fi p and 
fi<2 are varied in time. To this end, we adopt the linear 
stability analysis which has wide application in various 
nonlinear systems. 

To begin with, we expand the state vector as: 

tpi = V- + Srpt , 

where Stpi represents the deviation of the probability am- 
plitude in state \i) from the CPT solution ©. Adiabatic- 
ity is obeyed as long as these deviations remain small. 
Inserting this expansion into the dynamical Eqs. JSJ, we 
obtain the linearized equations: 



i— Sip = MSip + V , 

dt 



(7) 



where 



Sip 
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, M= | 







fi p < 

fi p v a ° 
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The last term at the right hand side of Q represents 
the "source" term arising from the temporal variation of 
Rabi frequencies (fi p and fid). 

The eigenfrequencies of the linearized equations (i.e., 
the eigenvalues of matrix M) can be easily found as 



wo = 0, 
1 

w ± = 77 



A±(A 2 + fi d fi^) 1/2 ' 



The corresponding eigenstates, i.e., the normal modes, 
are 



w 




w± 



:A4 




where 7Vo,± are normalization constants. Note that wo,± 
are orthogonal to each other. 

Since the frequencies u>o.± are an real, in the limit of 
fip,d - ► (i.e., the vanishing source terms), Stpi will not 
grow in time and hence the system initially prepared in 
a CPT state will stay in the instantaneous CPT state. 
This is the adiabatic theorem generalized to this non- 
linear system. We remark that for a general nonlinear 
system, complex eigenfrequencies of the linearized equa- 
tions may appear |12| . These systems are dynamically 
unstable such that, even in the absence of source terms, 
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the deviations may grow from intrinsic quantum fluctu- 
ations 0, [^} . Adiabaticity will therefore break down in 
the dynamical unstable regimes of a nonlinear system. 

The lack of complex eigenfrequencies in the current 
study guarantees that the system under consideration is 
always dynamically stable, hence the growth of devia- 
tions dipi can only be driven by the source terms. 

From the first sight, the presence of the zero mode 
wo seems to be most worrisome, since the fluctuation in 
this mode is energetically resonant with the CPT state. 
Fortunately, as we now show, this resonance does not lead 
to the growth of the zero mode. Writing the deviation 
vector as superpositions of normal modes 



The adiabatic condition now becomes 



Sip = 



E 

i=0, ± 



Equations J7J yield: 

ic a = uj a c a - i w£ * ■ 
In particular, we have 

ZC - -t W Q * = -iNo f - ^ V£ + $ 



Using the CPT solution given in Eqs. JO}, it is straight- 
forward to show that the right hand side of the above 
equation vanishes and hence the amplitude of the zero 
mode fluctuation does not change in time. Therefore we 
reach a very important conclusion: Although the zero 
mode fluctuation is resonant with the CPT state, it is 
not coupled by the dynamics. From now on, we can sim- 
ply ignore the zero mode and focus on the remaining two 
modes u>±. 

Assuming c±(0) = 0, the amplitudes c±{t) can be 
solved as 



c±(i) 



dt' e lu±{t '- 



(t')Vo(t'). (8) 



Defining the adiabaticity parameter as 
1 



r nl (t) = ~VM*)l 2 + M*)| 2 , 



(9) 



which represents the population in the fluctuations above 
the CPT solution, the adiabatic condition can now be 
defined quantitatively as 



r n i(t) < 1 



(10) 



When oj± is large, the exponential terms in JHJ is 
rapidly oscillating and the most significant contribution 
of the integral comes from t w /. We may change the 
time variable of -w±(t')*&o(t') in the integrand to t and 
take it out of the integral, which yields 



c±{t) = -wl(t)9 (t) 



l-e 



-A4 



10J± 



a, 



Qnl 
"off 



r„i(i) 



' + — I — — — <i 



L0_ 



a, 



^cff 



-iuj±t 



where we have neglected the oscillating term e 
evaluating r n \(t). 

For A = ; 0, after some algebra, the above inequality 
leads to 



r-ni = 



1 



< 1 



(11) 



10J± 



where \ = ^p/^d a s before. Comparing (|ll|l with we 
see that the adiabatic condition for the nonlinear system 
has a somewhat similar expression to that for the linear 
system, with, however, the following major difference: \ 2 
in the denominator at the left hand side of is replaced 
by \/l + £>X 2 i n CD- It follows that, at the later stage 
of STIRAP where \ ^ E the adiabatic condition for 
the nonlinear system is harder to fulfil than for its linear 
counterpart. 

We remark that for the linear system, we can use the 
same linearization procedure to derive the adiabatic con- 
dition. It is obvious that, in this case, the eigenfrequen- 
cies of the linearized equations are the same as the ones 
of the original equations. Eqs. J^J and will then 
reproduce the linear adiabatic condition J2J). 

To support the above analysis, we use the following nu- 
merical example. We choose two equal-amplitude equal- 
width Gaussian pulses for the pump and dump field cen- 
tered at t Pt d, respectively, 

where we have taken the pulse width as the units for 
time. A counter-intuitive pulse sequence |llj | requires 
that td < t p . At the initial time (t = 0), only the \a) 
state is populated, i.e., 

^ a (0) = 1, ^(0)=^m(0)=0. 

We evolve the state vector under the Hamiltonian Hu n 
and H\ n for the linear and nonlinear systems, respec- 
tively. 

Figure|21shows an example of the population dynamics. 
There are not much differences at the early stage for the 
linear and nonlinear cases. However, at the later stage of 
the STIRAP, the nonlinear system can no longer follow 
the CPT state while the linear system stays closely in the 
CPT state. In the end, we achieve an 80% population 
transfer in the former and near perfect transfer in the 
latter. 

Figure|21shows the corresponding time evolution of the 
adiabaticity parameter. At the early stage of STIRAP, 
the two parameters in the linear and nonlinear systems 
are almost identical. However, at the later stage, we have 

r n \ > Hin , 
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FIG. 2: (Color online) Population dynamics. Solid lines: pop- 
ulation in state \a) (\i/j a \ 2 ); dashed lines: population in state 
\g) (\4>g\ 2 for the linear case and 2|^/> g | 2 for the nonlinear case). 
The parameters used here are A = 0, flo ~ 5, td ~ 3 and 
tp = 3.8. 




time 

FIG. 3: (Color online) Adiabaticity parameter as defined in 
© and llll l for the linear and nonlinear systems, respectively. 
Same parameters as in Fig. 



in perfect agreement with our analysis, which also ex- 



plains the reduced population transfer efficiency for the 
nonlinear system and the population dynamics as demon- 
strated in Fig. 

In conclusion, we have shown that, for nonlinear sys- 
tems, the textbook method used in deriving the adiabatic 
condition fails to provide the correct answer, and the 
proper way to derive the adiabatic condition is through 
a linearization procedure. This is demonstrated by con- 
sidering the coherent population transfer process via STI- 
RAP in both a linear and a nonlinear three-level A- 
system, which are indeed governed by different adiabatic 
conditions. As a side, our analytical result of the non- 
linear adiabatic condition will become very useful in 
determining the required laser parameters for photoasso- 
ciation via STIRAP. 

As we have mentioned, in order to keep the mathemat- 
ics as simple as possible, we have neglected here the non- 
linear collisions between particles which are important 
in the ultracold quantum degenerate atomic/molecular 
systems. The inclusion of the collisions will introduce 
dynamically unstable regimes for the CPT states [T^| . 
With collisions, the linearization procedure is equivalent 
to the Bogoliubov treatment of weakly interacting Bosc 
condensates. However, the orthonormality of linearized 
eigenstates is replaced by bi- orthonormality of the Bo- 
goliubov normal modes, and the zero-frequency mode 
corresponds to the Goldstone modes related to the spon- 
taneously broken symmetries. These changes lead to con- 
siderably more complicated algebra. The adiabatic con- 
dition for collisional systems will be studied elsewhere. 
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